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SUMMARY 

J 

This  Memorandum  describes  work  done  to  improve  the  measure  of  agreement 
between  a  target  one-third-octave  spectrum  and  its  reproduction  via  some  remote 
sound  source.  Digital  filters  create  limited  bandwidth  data  which  are  modified 
and  finally  recombined  to  produce  a  new  time  history.  Studies  of  certain  fixed 
and  rotary  wing  aircraft  have  shown  mean  errors  significantly  reduced  by  this 
technique . 


Copyright 

© 


Controller  HMSO  London 
1987 


2 


LIST  OF  CONTESTS 

1  INTRODUCTION 

2  METHOD 

2.1  Principle  of  method 

2.2  Analysis  by  digital  filter 

2.2.1  General 

2.2.2  Performance 

2.3  Synthesis 

2.4  Co’Tfe^tion  file 

2.5  Postcript 

3  CONCLUSIONS 
Acknowledgment 

Appendix  A  General  theory 

Appendix  B  Derivation  of  band-pass  filter  used 

Appendix  C  Iteration  process  (system-specific  details) 

References 

Illustrations 

Appendix  Illustrations 

Report  documentation  page 


Page 

3 

4 

4 

5 

5 

6 

6 

7 

7 

8 
8 
9 

14 

19 

20 

Figures  1-4 
Figures  A1-A8 
inside  back  cover 


TM  FS(F)  679 


3 


1  INTRODUCTION 

A  requirement  in  many  ground-based  experiments  is  the  ability  to  establish 
a  determinate  sound  pressure  field  in  a  chamber  that  can  accommodate  human  sub¬ 
jects.  Frequently,  this  field  is  required  to  approximate  to  that  existing  in 
some  distant  environnent,  such  as  an  aircraft  cockpit  in  flight. 

A  tape-recording  of  the  time-history  of,  say,  the  cockpit  noise  can  be 
obtained  and  may  well  represent  to  the  ground-based  observer  the  only  available 
information  on  the  frequency  content  of  the  source  field.  We  may  conveniently 
refer  to  this  as  the  target  field  in  both  time  and  frequency  domains. 

However,  on  replay  in  the  ground  chamber,  that  tape-recording  will  not  in 
general  produce  an  exact  duplication  of  the  target  field  for  many  reasons:  the 
multiple  mechanical  sources  of  the  original  noise  and  the  acoustics  of  the  cock¬ 
pit  cannot  be  precisely  reproduced  in  a  laboratory.  Thus,  frequency  analysis 
equipment  set  up  in  the  chamber  during  the  time  that  some  sound-reproduction 
equipment  is  replaying  the  target  time  history,  will  not  display  the  target 
spectrum  as  existed  in  the  aircraft,  but  rather  an  observed  spectrum  that  differs 
from  it  by  an  error  spectrum. 

A  typical  third-octave  error  spectrum  is  shown  in  Fig  la  for  a  modern  jet 
aircraft  where  the  target  field  is  that  of  whitish  noise  over  the  audible  range 
of  the  spectrum;  it  is  seen  that  deviations  of  ±10  dB  exist, 

A  method  has  been  developed  whereby  the  error  spectrum  may  be  reduced  con¬ 
siderably  as  seen  in  Fig  lb.  The  method  is  lengthy,  even  with  the  use  of  an 
Array  Processor,  but  it  has  proved  practicable  in  dealing  with  typical  time 
histories  of  aircraft  noise  where  the  duration  of  some  quasi-stable  state  of 
interest  does  not  exceed  30  seconds  or  so. 

Fig  IcSd  are  a  comparable  pair  of  error  spectra  for  a  noise  source  with 
significant  deterministic  frequencies  present  -  they  derive  from  a  helicopter  at 
the  hover.  Again,  a  considerable  improvement  is  noted,  this  time  extending  down 
to  the  low  harmonics  of  the  blade-passing  frequencies. 

It  has  been  observed  that  the  revised  aircraft  noise  markedly  retains  the 
character  of  the  original  noise,  as  opposed  to  that  of  totally  synthetic  random 
noise . 
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2  METHOD 

2.1  Principle  of  method 

The  principle  of  the  method  is  as  follows.  An  observed  time  history  is 
obtained  by  means  of  a  microphone  in  the  chamber,  and  this  is  analysed  in  real 
time  by  a  B/K  frequency  Analyser  Type  2131.  The  resultant  data  file,  covering 
third-octave  bands  from,  say,  10  Hz  to  8  kHz  is  subtracted  from  a  corresponding 
data  file  obtained  by  identical  hardware-analysis  of  the  target  time  history. 

This  difference  file  forms  the  basis  for  the  (first)  error-correction  file. 

Then  follow  two  processes,  broadly  effecting  analysis  of  the  error,  then 
synthesis  of  a  compensated  time  history  that  will  replace  the  original  aircraft 
recording. 

In  the  analysis  module,  the  target  time  history  (stored  on  disc)  is  applied 
to  a  recursive  digital  filter  and  decomposed  into  its  'component  band-files'  - 
this  is  a  Fourier-like  process  in  the  time  domain  producing  one-third  octave-wide 
files. 

The  component  band-files  (which  may  be  numerous)  are  stored  on  disc  pending 
synthesis  by  a  complementary  module.  This  produces  a  modified  time  history  com¬ 
prising  the  sum  of  the  component  band-files  weighted  by  coefficients  that  are 
derived  from  the  values  in  the  error  correction  file. 

At  the  end  of  the  process,  a  new  time  history  file  will  exist  that  contains 
all  the  constituent  frequency  components  of  the  target  file,  but  these  being 
differentially  weighted,  a  different  noise  field  will  be  produced  on  playing  the 
new  file  into  the  chamber. 

If  the  resultant  new  (ie  second)  error  spectrum  is  not  sufficiently  small 
it  may  be  used  to  produce  a  further  weighted  sum  of  component  band-files  via  a 
repeat  application  of  the  synthesis  module,  and  so  on.  However,  small  differ¬ 
ences  in  chamber  topography,  amplifier  gains  etc  impose  an  early  limit  to  the 
number  of  iterations . 

Existing  procedures  in  the  laboratory  enable  data  to  be  acquired  at  a  high 
sampling  rate  (32  kHz)  and  written  out  co  300  mB  disc  storage.  Other  procedures 
allow  disc  data  to  be  put  out  at  this  rate  into  a  high-power  sound  reproduction 
system.  Fig  2  -shows  in  principle  the  mechanisation  of  the  whole  package. 
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2.2  Analysis  by  digital  filter 


2.2.1  General 


The  process  of  filtering  the  target  waveform  is  thus  fundamental  co  the 
method.  The  filter  used  is  effectively  a  Buttervorth  (maximally  flat)  4th-order 
bandpass  filter  with  zero  phase  shift,  implemented  on  the  Array  Processor  (AP) 
and  windowed  so  that  no  discontinuities  are  introduced  although  only  a  small 
portion  of  the  whole  time  history  can  be  accommodated  for  processing  at  any  one 
time  in  the  AP. 

Rather  than  merely  define  the  filter  in  terms  of  a  difference  equation,  or 
alternatively  as  a  z-plane  transfer  function,  it  has  been  thought  worthwhile  to 
outline  the  derivation  of  these  expressions  in  Appendix  B,  which  itself  seemed  to 
necessitate  Appendix  A,  where  an  outline  is  given  of  some  relevant  features  of 
s-plane  and  z-plane  analyses. 

The  solution  to  a  difference  equation  having  the  form  of  equation  (B-9) 
was  conveniently  already  existing  as  a  module  in  the  APMATH  library  of  the  Array 
Processor.  As  shown  in  Appendix  B,  the  band-pass  filter  derives  from  a  first- 
order  lag  transfer  function. 

After  passing  forwards  through  the  filter,  the  data  is  time-reversed  for 
a  second  pass  which  cancels  the  phase  shift  introduced  on  the  first  pass',  while 
providing  steeper  attenuation  characteristics  in  the  stop  bands.  The  forward/ 
reverse  passes  are  then  repeated  for  further  selectivitiy . 

As  the  object  of  the  process  is  finally  to  synthesize  a  revised  time  history 
without  added  discontinuities,  only  the  central  portions  of  each  data  set  being 
filtered  are  in  fact  retained.  To  elaborate:  storage  in  the  AF  allows  32K  data 
points  (te  representing  1  second  of  time  history)  to  be  processed  at  one  time, 
allowing  for  the  necessary  location  shifts  in  memory  during  manipulation  of  the 
data.  The  32K  data  points  are  first  windowed  by  cosine  tapering  to  zero  of  the 
first  and  last  4K  points,  and  then  only  the  central  16K  data  points  are  retained 
as  a  {  second  output  component  band-file.  Thus  n  seconds  of  target  file  time 
history  produce  (2n  -  1)  component  band-files  each  J  second  long. 

As  the  data  passes  repeatedly  through  it,  the  filter  has  its  third-octave 
pass  band  centred  sequentially  on  the  standard  ISO  centre  frequencies,  to  be  com¬ 
patible  with  the  B/K  2131  for  hardware  analysis:  hence  for  decomposition  into  m 
1/3-octave  bands,  there  will  be  produced  m(2n  -  1)  J-second  files,  stored  on 
disc.  Clearly,  this  must  be  a  time-consuming  process. 
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2.2.2  Performance 

The  attenuation/frequency  characteristic  of  the  filter  will  now  be  discus¬ 
sed  briefly.  In  fact,  it  is  preferable  to  speak  of  the  aforesaid  characteristic 
as  pertaining  to  the  process,  rather  than  to  the  filter,  because  as  has  been 
seen,  the  output  of  the  process  is  a  summed  function  of  individual  filter  out¬ 
puts.  Thus  instead  of  exciting  a  single  filter  by  stepping  through  a  succession 
of  sinusoids  and  plotting  the  output  attenuation  in  one  band,  a  more  valid  result 
is  obtained  by  plotting  the  simultaneous  output  of  all  filters,  for  one  forcing 
sinusoid. 

Accordingly,  Fig  3a  shows  the  response  to  a  single  (ta  FWD  +  REV)  pass 
through  the  process  for  a  1  kHz  input.  It  is  seen  that  the  skirts  of  the  filter 
response  fall  off  at  about  -12  dB/octave,  but  the  initial  rate  is  about  -33dB/ 
octave.  The  first  figure  is  compatible  with  the  -6  dB/octave  characteristic 
expected  from  each  pass  through  the  filter  (ie  FWD),  remembering  that  the  second 
(REV)  pass  multiplies  the  attenuations  and  produces  a  further  -6  dB  in  the  final 
result.  It  is  not  known  why  the  initial  rate  of  attenuation  is  as  high  as  has 
been  observed  above.  Fig  3b  shows  the  overall  characteristic  vhen  a  double  pass 
(ie  2  x  (FWD  +  REV))  is  used.  As  expected,  the  logarithmic  slopes  are  doubled. 

2.3  Synthesis 

Conceptually,  this  is  simply  the  concatenation  of  all  the  J  second  compon¬ 
ent  band-files  already  stored  on  disc,  and  a  summation  over  all  bands  to  produce 
the  final  output  file.  In  fact,  the  data  is  handled  again  by  the  AP,  in 
1-second  slices  -  the  required  weighting  factors  are  applied  and  a  sequence  of 
multiply-and-add  operations  builds  up  in  the  AP  a  composite  J-second  file  with 
weighted  contributions  from  all  its  component  band-files.  This  1-second  output 
file  is  then  sent  back  to  be  concatenated  on  disc  with  its  issnediate  predecessor 
in  time.  Thus  the  length  of  the  final  output  file  is  (2n  -  l)/2  seconds,  ie  half 
a  second  shorter  than  the  original  target  time  history. 

A  special  case  of  the  output  file  is  that  where  the  weightings  are  unity  - 
when  a  file  whose  characteristics  match  that  of  the  original  target  file  should 
be  produced.  Fig  4a  shows  a  portion  of  a  time  history  of  a  helicopter  noise 
field,  and  Fig  4b4c  the  reconstruction  by  such  unit  weighting  for  single  and 
double  complete  passes  respectively.  It  should  be  noted  that  because  of 
undesired  but  unavoidable  attenuation  in  the  pass  band  of  the  filter,  the  fre¬ 
quency  components  present  in  the  target  file  (which  are  unlikely  to  lie  at  exact 
band-centre)  have  amplitude  errors  which  result  in  the  reconstructed  waveforms 


not  precisely  imaging  the  original.  For  this  reason,  the  waveforms  of  Fig  4  do 
not  have  identical  y-axis  scales  and  these  are  omitted  from  the  plots,  which  are 
intended  to  be  of  a  qualitative  nature.  It  has  in  fact  been  found  to  improve  the 
iteration  process  if  the  first  correction  file  produced,  to  be  used  in  producing 
the  first  iterated  synthetic  output  file,  is  actually  derived  via  this  'unity- 
reconstructed'  file  and  not  the  true  target  file  (which  of  course  still  provides 
the  ultimate  error  spectrum  to  be  minimised)  . 

2.4  Correction  file 

It  should  be  mentioned  that  before  the  laboratory  computer  system  was  modi¬ 
fied  to  include  on-line  hardware  spectral  measurement  and  data  storage  via  the 
B/K  2131,  considerable  effort  was  put  into  developing  off-line  routines  that 
would  accept  narrow-band  (1  Hz)  spectral  data  via  FFT  programs  in  order  to  pro¬ 
duce  1/3-octave  data  to  be  used  in  the  derivation  of  correction  files.  The  pro¬ 
cess  clearly  adds  even  more  time  to  that  already  needed,  and  it  was  necessary  to 
model  the  I /3-octave  band  filters  in  the  B/K  equipment  fairly  closely,  because 
in  the  final  analysis,  a  verdict  on  the  goodness  of  the  compensation  was  still 
made  via  on-line  measurement  using  the  B/K  2131.  As  stated,  the  software  1/3- 
octave  routine  proved  unnecessary  in  the  end. 

Other  approaches  to  the  optimisation  problem  have  been  tried  with  little 
success,  that  is,  by  using  more  complicated  algorithms  to  derive  a  correction 
file,  rather  than  merely  getting  weighting  factors  from  the  error  file.  The 
difficulty  is  that  the  decomposition  process  of  section  2.2.1  is  obviously  far 
from  perfect  in  both  pass  and  stop  bands,  and  it  was  not  found  possible  to  reduce 
errors  to  acceptably  small  limits  by  making  use  of  the  estimated  imperfections  in 
the  filtering  process.  As  already  stated,  repeated  derivation  of  error  files  and 
their  application  as  correction  coefficients  has  been  found  to  be  the  only  way  to 
achieve  the  improvements  noted  {eg  Fig  J). 

2.5  Poster ipt 

A  revised  modus  operandi  has  been  developed  since  the  main  text  was 
written,  and  brief  notes  are  given  here. 

It  has  always  been  appreciated  that  there  was  a  problem  in  dealing  with 
che  lower  frequency  bands  that  are  of  especial  interest  in  handling  helicopter 
noise.  For  example,  at  a  centre  frequency  of  20  Hz  the  1/3-octave  bandwidth  is 
about  ±2.5  Hz:  the  effective  truncation  of  analysis  time  to  1  second  means  that 
even  after  the  retention  of  only  the  centre  half  second  and  use  of  tapering,  the 
transient  (impulse)  response  of  the  filter  to  the  too-short  section  of  time  his¬ 
tory  can  be  recognised. 
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Ac  the  same  time,  difficulties  were  being  met  in  trying  to  implement  alter¬ 
native  filters  of  higher  order.  Round-off  computational  errors  have  tradition¬ 
ally  restricted  the  filter  order  to  one  or  two,  but  here  it  was  found  thac  as  the 
centre  frequency  and  bandwidth  decreased,  the  second-order  solution  became 
increasingly  erroneous  and  finally  unstable  for  centre  frequencies  typically 
around  50  Hz.  Double  precision  arittxnecic  had  already  been  used  as  far  as  poss¬ 
ible  when  it  was  realised  that  the  use  of  the  AP  in  the  analysis  module  was  the 
cause  of  the  trouble.  Although  the  AP  word-length  is  greater  than  the  single 
precision  word- length  on  the  PDP  11/84,  it  is  much  shorter  than  the  double  pre¬ 
cision  word-length;  thus  in  effect  accuracy  was  being  sacrificed  to  speed. 

Accordingly,  the  program  has  been  reconfigured  for  double  precision  working 
completely  on  the  PDP  11/84,  with  the  very  considerable  advantage  that  there  is 
no  longer  any  need  to  analyse  in  steps  of  1  second  -  it  becomes  a  continuous  pro¬ 
cess  over  the  whole  duration  of  the  input  time  history.  The  AP  is  only  (and  very 
advantageously)  used  in  the  synthesis  module  where,  as  before,  it  functions  as  an 
accumulator,  but  this  time  the  process  is  considerably  speeded  up,  because  the 
amount  of  file  manipulation  on  disc  is  reduced.  The  overall  time  for  analysis 
and  synthesis  is  of  the  same  order  as  originally. 

•3  CONCLUSIONS 

A  digital  method  has  been  demonstrated  for  analysing  a  target  signal  field 
and  synthesizing  another  one  related  to  it.  The  method  is  simplistic,  lengthy 
and  approximate,  but  has  been  successfully  mechanised.  Significant  improvements 
have  been  observed  in  the  reproduction  of  aircraft  cabin  noise  fields. 

It  is  expected  that  further  work  will  be  carried  out  using  digital  band¬ 
pass  filters  with  different  frequency  characteristics  with  a  view  to  effecting 
optimum  performance. 
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Appendix  A 


GENERAL  THEORY 


The  s-plane  (complex  frequency)  has  long  been  deemed  to  be  of  more  utility 
than  the  (real)  time  plane  as  an  aid  to  the  design  and  analysis  of  Continuous 
Control  Systems  (in  the  widest  interpretation  of  their  definition) . 

The  Laplace  Transform  F(s)  of  a  time  function  f(t)  allows  direct  access 
to  the  convenient  representation  of  modal  characteristics  and  system  responses  on 
the  s-plane.  The  alternative  manipulation  of  f(t)  using  convolution  is  in 
general  more  difficult  and  less  informative. 

When  sampled-data  He  digital)  systems  are  considered,  a  continuous  time- 
signal  f(t)  ,  sampled  every  T  seconds,  is  represented  by  f*(t).  The  latter 
has  a  Laplace  Transform  F*(s)  that  contains  non-algehraic  factors  exp(-nsT) 
for  n  *  0,1,2,,,.  and  is  not  very  tractable  in  this  form. 


A  convenient  transformation  of 


exp(sT) 


(A- 1  a) 


Chen  enables  the  Laplace  Transform  to  be  written  in  terms  of  z  ,  when  it  is 
known  as  the  z-Transform.  In  this  form  the  awkward  factors  referred  to  above 

become  _ 

z  n  for  n  »  0,1,2,... 

and  they  can  be  interpreted  as  Delay  operators  in  the  transformed  time  series, 

F  (z)  -  e(0)  +  e(T)z"’  +  e(2T)z'2  +  e(3T)z“3  +  ...  . 

This  method  of  describing  the  sampled  time  series  f*(t)  in  terms  of  the 
delay  operator  z  is  very  powerful,  enabling  transfer  functions  of  discrete 
processes  to  be  expressed  as  rational  functions,  and  difference  equations  to  be 
derived  from  them,  giving  the  current  output  value  y(i)  as  a  sum  of  past  output 
values  y(j)  for  j  »  i  -  1  ,  i  -  2,...  and  present  and  past  input  values  x(k) 
for  k  »  i  ,  i-l,i-2,....  This  will  be  done  in  Appendix  B  for  the 

filter  under  discussion. 

However,  it  will  be  realised  that  the  above  definition  of  z  in 
equation  (A- la)  implies  the  inverse  relation 


logg(z) 


(A- lb) 


and  if  the  digital  equivalent  H(z)  to  a  continuous  transfer  function  H(s)  is 
sought,  then  some  approximation  to  the  logarithm  must  be  used.  A  well-known 
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approximation  leads  to  equation  (A-lc)  in  the  next  section,  but  first  there 
follows  a  brief  introduction  to  some  graphical  techniques. 

A. I  Mappings 

The  fundamental  feature  of  sampled  data  is  that  unless  a  continuous  time 
series  is  band-limited  to  below  the  Nyquist  frequency  f^  «  f^/2  ”  d/2T  Hz), 
then  no  longer  can  the  process  be  uniquely  inverted  to  reproduce  the  original 
waveform.  This  mechanism  (aliasing)  occurs  because  of  the  duality  existing 
through  the  Fourier  Transform  between  the  processes  of  multiplication  and  convo¬ 
lution  in  the  time  and  frequency  planes. 

Thus,  because  a  time  series  is  in  effect  a  sequence  of  products  of  unit 
delta  function  pulses  and  a  continuous  time  history  (Fig  AI)  then  in  the  fre¬ 
quency  plane  the  resultant  convolution  is  seen  to  Be  periodic  with  repetitions 
of  the  spectrum  at  intervals  of  f  (see  Ref  3) . 

Hence  in  Fig  Al  the  part  of  the  spectrum  containing  genuine  information  is 
the  Base-band  from  -fs/2  to  *fg/2  and  clearly  this  will  be  corrupted  by 
spurious  folding-down  of  frequencies  from  higher  bands  unless  the  above-mentioned 
band  limiting  is  assured. 

Looking  at  the  (discrete)  s-plane,  it  is  seen  that  the  area  of  interest  is 
confined  to  the  strip  between  s  *  jw^/2  and  s  ”  -jwg/2  .  Any  poles  or  zeros 
appearing  in  this  region  (Fig  A2)  will,  by  virtue  of  the  periodicity  noted  in 
Fig  Al ,  also  appear  in  all  other  complementary  strips. 

As  the  imaginary  axis  forms  the  boundary  between  stable  and  unstable  modes 
of  system  behaviour,  design  attention  in  the  discrete  s-plane  is  focussed  on  the 
so-called  Primary  Strip  in  Fig  A2. 

Recalling  the  transformation  z  »  exp(sT)  it  can  be  shown  readily  by  sub¬ 
stituting  s  *  o  +  jiu  that  the  primary  scrip  maps  into  a  circle  on  the  z-plane, 
see  Fig  A3a.  The  radius  of  the  circle  is  unity  and  the  mappings  of  a  single  pole 
and  a  pair  of  complex  poles  is  also  shown  in  Fig  A3a  for  clarity.  It  is  seen 
that  poles  outside  the  unit  circle  will  result  in  unstable  modes. 

As  stated  in  the  preceding  section,  much  design  experience  in  the  field  of 
continuous  (analogue)  filters  can  be  made  use  of  by  an  approximation  to  the 
z-transform,  which  is  equivalent  to  mapping  the  unit  circle  on  the  z-plane  to  a 
new  plane.  Different  authors  have  used  different  forms  of  this  transformation 
with  the  result  that  their  terminologies  and  illustrative  diagrams  can  cause 
confusion. 
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Essentially,  the  transformation  is  the  bilinear  relation 


1 


*  Jv,  .  1SL±  L  .  'll!!)  . 

1  1  Z  +  1  ^  W|  -  1  J 


(A-lc) 


This  mapping  projects  the  unit  circle  on  the  z-plane  into  the  entire  LH 
Wj-plane  (Fig  A3b) .  It  can  be  seen  that  traversing  the  unit  z-circle  from  (I) 
to  (2)  corresponds  to  the  whole  axis  of  positive  frequencies  on  the  new  plane. 
If,  therefore,  a  mapping  is  required  between  the  entire  LH  continuous  Wj-plane 
and  the  primary  strip  of  the  discrete  s-plane  (Fig  A3a) ,  via  the  z-plane,  then 
the  relation  between  the  corresponding  frequencies  must  be  highly  nonlinear,  or 
warped . 

Three  variants  of  the  bilinear  transformation  are  commonly  found  and  will 
be  illustrated  briefly. 

(1)  This  is  the  basic  form  already  introduced  ahove  in  equation  A-lc). 
From  this  and  the  definition  z  *  exp(sT)  we  can  derive  the  relation 


'1 


tanh(sT/2) 


Then  by  setting  w  =  jv.  and  s  *  jw  we  get  the  warping  equation 


tan(u)T/ 2) 


(A- 2a) 


(A- 2b) 


which  links  a  given  frequency  in  the  discrete  s-plane  to  some  other  frequency  on 
the  continuous  Wj-plane. 

This  important  curve  is  shown  in  Fig  A4a  (see  Ref  4;. 

(2)  Correct  dimensionalising  and  normalisation  with  respect  to  the  sampling 
frequency  is  introduced  by  the  factored  mapping 

_  2(z  -  1) 


T(z  +  1) 


(A-3) 


which  leads  to  the  corresponding  expression 


and  frequency  mapping 


Y  tanh(sT/2) 


’2  -  y  tan(“T/2) 


(A-4a) 


(A-4b) 


This  is  shown  in  Fig  A4b,  drawn  to  illustrate  the  relation  between  frequency 
ratios,  as  can  be  found  in  the  literature^.  Also  shown  is  the  line  defining 
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This  is  equivalent  to  using  the  bilinear  transform  with  immediate  substitution 
of  continuous-plane  parameters,  and  is  seen  to  be  a  reasonable  approximation  to 
the  tangent  function  up  to,  say,  wT  •  1  .  Beyond  this,  frequency  warping 
becomes  significant  and  the  modified  continuous-plane  frequencies  should  be  used. 

It  is  important  to  realize  that  equation  (A-3)  is  an  exact  relation,  having 
certain  properties,  desirable  or  otherwise,  but  is  only  an  approximation  to  the 
exact  inverse  of  z  ■  exp(sT)  viz: 

*  “  f  ioggCz)  • 


The  right-hand-side  of  equation  (A-3)  is  in  fact  the  first  term  in  an  approxima¬ 
tion  to  log  (z)  as  was  mentioned  in  equation  (A-l),  or,  equivalently,  derived 
by  similarly  truncating  exp(sT)  ,  ie  exp(sT/2)/exp(-sT/2)  . 

Yet  another  way  of  visualising  equation  (A-3)  as  being  an  approximation  is 
by  considering  the  operation  (1/s)  -  integration  in  the  complex  plane  -  and  its 
counterpart  in  the  time  domain. 


Fig  A5  shows  that  when 

y  »  f  x dt 

is  evaluated  by  trapezoidal  approximation,  the  current  estimate  of  y  {ie  y  ) 


is  given  by 


Vi 


+  2  (xn  +  Vl> 


-1 


or,  recognizing  z  as  the  unit  delay  operator, 

Y(z)  -  Y(z)z_l  +  |  X(z)z  +  |  X(z)z"' 
which  gives  the  z-transfer  function  as 


Y(z) 


xTzT 


T(1  +  z~') 
2(1  -  z"') 


(A-5) 


This  expression,  describing  the  operator  (1/s)  ,  is  seen  to  be  correctly  the 
reciprocal  of  the  relation  already  given  in  equation  (A-3)  for  the  operator  (s)  . 

If  the  filter  is  being  designed  on  the  z-plane  with,  say,  cut-off  frequency 

u)  .  ,  the  substitution  of  w-  by  relation  equation  (A-3)  will  produce  a 

dc  z 

continuous-plane  approximation  where  the  cut-off  frequency  of  the  equivalent  con¬ 
tinuous  filter  is  at  a  frequency  id  .  Thus  we  could  have,  for  example,  the 

ac 

mapping  as  sketched  in  Fig  AAc,  enabling  u  to  be  read  off  against  a  known 
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(3)  A  further  development  of  (1)  is  the  transform 


where  K 


where  Sl^  is  a  (cut-off)  frequency  on  the  continuous  plane  for  which  there  is 
some  a  priori  reason  for  mapping  to  a  specific  similar  point  on  the  discrete 
plane  -  generally  because  tables  may  exist  for  the  design  of  a  particular  analogue 
filter.  For  generality  we  retain  the  parameter  !5cq  but  usually  it  may  be 
equated  to  unity. 

Then,  the  constant  K  (which  does  not  affect  the  shape  of  the  mapping) 
ensures  a  convenient  system  of  scales  for  the  graphs.  Thus 


Im(di^) 


tan  u  T7 
c 


as  shown  in  Fig  A4d  (see  Ref  6) . 
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DERIVATION  OF  BAND-PASS  FILTER  USED 

B. 1  General 

IC  has  been  shown  above  thae  methods  exist  whereby  the  continuous  s-plane 
may  be  mapped  to  the  discrete  s-plane.  In  particular,  we  must  now  derive  a 
z-transfer  function  for  a  band-pass  digital  filter,  so  that  it  may  be  implemented 
on  the  AP. 

The  steps  in  the  process  will  be: 

(1)  to  effect  a  supping  on  the  continuous  plane  between  a  prototype  low- 
pass  filter  and  a  general  band-pass  filter. 

(2)  To  extend  this  mapping  into  the  z-plane  by  obtaining  an  approximate 
relation  between  s  and  z  that  will  map  the  prototype  low-pass  filter  into  a 
digital  band-pass  filter,  without  determining  its  numerical  parameters. 

(3)  To  obtain  such  parameters  by  substitution  of  z  »  exp(sT)  into  the 
relation  found  in  (1)  above. 

(4)  To  obtain  the  z-transfer  function  by  substituting  these  numerical 
values  in  (2)  above  and  thence  write  down  the  difference  equation. 

B. 2  Details 

(1)  Guillemin  has  shown^  that  other  more  compli.ated  continuous-plane 
filters  may  be  mapped  by  frequency  transf ormations  from  the  low-pass  form  which 
is  regarded  as  a  prototype.  (This  preparatory  mapping  solely  in  the  s-plane  is 
not  to  be  confused  with  s-  to  z-plar.e  mappings  considered  above.  ) 

Take  the  prototype  , 

pa)  '  —  (175  1  *  (B-,) 

CO 

where  X  *  a  ♦  jw 
B  rJ 

anc*  ^co  *■*  frequency. 

(This  nomenclature  is  chosen  arbitrarily  to  allow  the  familiar  complex  variable 
s  *  a  *  ju  to  be  used  as  the  output  variable.) 

We  require  the  new  variable  s  to  be  such  that  given  P(X)  as  in 
equation  (B-l)  above,  then  P(s)  will  describe  a  band-pass  function.  Values  of 
X  in  equation  (B-l)  are  next  mapped,  noting  that  in  the  X-plane 
X  ■  0  gives  the  centre  of  the  (low)  pass  band, 

X  -  •>  gives  the  high-frequency  asymptotic  behaviour. 


L 
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(B-3a) 
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[In  passing,  ve  note  that  by  substituting  equation  (B-3a)  into  equation  (B-l)  we 
obtain  the  s-plane  transfer  function  of  the  band-pass  filter  in  terms  of 
and  ,  viz 


8(“ u  •  “i> 


s  +  s(uiu  -  u>4)  ♦ 


( B—  3b  > 


Fig  A7  shows  the  familiar  shape  of  the  transfer  function  of  equation  (B-3b); 
however,  it  is  equation  (B-3a) ,  the  frequency  transformation  function,  that  will 
be  used  in  the  next  section.] 

(2)  The  mapping  of  the  prototype  low-pass  filter  (X-plane)  onto  the  z-plane 
is  via  the  s-plane  using  equation  (B-3a) .  Thus,  using  the  bilinear  transform 
by  putting  s  «  k(z- !)/(«  ■*•!)  in  equation  B-3a)  ve  get 


(z 


2Az~*  +  1 


-2 


(B-4) 


where 

and 


(k  ♦  a)  > 

n  . . .  B 

co  kin  , 

,2-.  4 

t— 

k  ♦  <ii 


1  -  z 


,  using  the  nomenclature  u> 


u  111. 
u  4 


(Note  that  k  is  still  undetermined,  and,  in  fact,  as  the  coefficients  A  and 
B  are  those  actually  required,  k  will  never  be  found  explicitly.) 


This  expression  for 
to  give 


Y(z)  a 
XT zT 


X  is  now  substituted  into  the  prototype  equation  (B-l) 

Will 


-2 


FM 


-2  -1 
z  -  2Az  ♦ 


(i) 


(B-5) 


+  1 


(3)  To  determine  the  constants  A  and  B  in  section  (2)  above,  we  next 
substitute  for  powers  of  z  in  equation  (B-4)  remembering  that  for  points  on 
the  unit  circle  on  the  z-plane  we  have 

z  *  *  exp(j(uT)  »  cos(iiiT)  -  j  sin(wT) 

Some  algebra  then  reduces  the  right-hand  side  of  equation  (B-4)  to 


jB(A  -  cos(uiT) 
sin(uT) 


(B-6) 


This  maps  to  the  general  continuous  frequency  u>  «  Im(A)  ,  and  we  proceed  with 

D 

the  mapping  of  other  frequencies  thus: 
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we  gee 


\  *  0  (mapping  to  band-pass  centre-frequency) 


A  «  cos(uT)  , 


“00  "  ?  C0S~'A  * 

To  use  the  mapping  at  the  -3  dB  points,  referring  back  to  Fig  A6  we  put 
Im(\)  »  uig  ■  «(1  and  *^co  successively,  thus: 

B(A  -  cos(w.T)) 

-si  -  - - 

co  sinfu.T) 

and 

B(A  -  cos(u  T)) 

+  SI  —  1  ■■  i  ■  v  —  r;v.— ■  ■ 1  -  ■ 

co  sinCw^T) 

From  these  equations  we  get  the  following  relations  - 
by  division 

cos  [U  *  »  )  |] 

A  .  - p - L_|i  ,  (E 

cos  [(u,u  -  Ujl) 

By  subtraction 


tanL(uu "  V  ij 

Thus  the  coefficients  in  equation  (B-5)  may  now  be  calculated  to  enable  the  final 
difference  equation  to  be  written  down.  ^co  is  arbitrary  and  for  convenience 

is  set  equal  co  1  rad/s.) 

(4)  Using  the  delay  property  (shifting  theorem)  of  the  bilinear  transform 
on  equation  (B-5)  gives  the  time  domain  solution 

(l  ♦  i)y.(t)  -  2Ay._,(t)  ♦  (l  -  “  fi)Xi(t)  ’  (l)  Xi-2(t) 


yi(c)  *  ttt-bt  *i<c>  -  o4tt  xi-2(t)  +ir^BTyi-i(t) 

+  irTTy  yi-2(t>  •  (B-9) 

This  is  the  equation  which  is  implemented  on  the  Array  Processor  after  the 
coefficients  have  been  inserted  from  equations  (B-7)  and  (B-8) . 
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It  is  to  Be  noted  that  the  relation  of  equations  (B-4)  and  (B-6)  is  shown 
plotted  in  Fig  A.8.  This  is  the  warped  relation  between  frequencies  on  the  con¬ 
tinuous  plane  for  the  prototype  lowpass  filter  and  those  on  the  discrete  plane 
for  the  bandpass  filter.  This  agrees  veil  with  Fig  3a  which  was  shown  for  a 
forward  pass  followed  By  a  reverse  pass. 


ITERATION  PROCESS  (SYSTEM- SPECIFIC  DETAILS) 

For  local  convenience,  there  follow  working  notes  on  the  procedures  that 
must  be  followed  to  go  through  an  iteration.  It  is  assumed  that  a  good  quality 
tape  recording  exists  (preferably  in  PCM  format  to  retain  maximum  S/N  ratio, 
together  with  the  system  hardware  as  outlined  in  Fig  2) . 

A  data  analysis  package  ’DATS'  (Prosig  Ltd)  is  present  on  the  system  and 
its  modules  and  facilities  have  been  used  extensively  in  the  work.  The  procedure 
'TESTX4'  that  puts  out  disc  data  in  real  time  to  the  analogue  power  amplifiers 
is  part  of  a  suite  of  programs  developed  by  Signal  Computing  Ltd. 

(1)  Acquire  onto  300mB  disc  a  suitable  recording  of  target  noise  using  PCM 
recorder  and  DATS/ACQUIR  . 

(2)  Use  DATS/ARITH  to  scale  the  numbers  in  (1)  to  a  level  suitable  for  sub¬ 
sequent  output  via  the  A/D  convertor  to  the  power  amplifiers. 

(3)  Obtain  the  1/3  oct.  spectrum  of  the  0/P  file  in  (2)  using  B/K  Type  2131 
Analyser.  Use  /AC2131  to  acquire  this  data  onto  disc  file. 

(4)  Run  LASSEX  (a  DATS-type  JOB  that  calls  the  digital  filter  module  /LASSTP  ) 
on  the  target  noise  time  history  obtained  in  (2). 

(5)  Use  /XASCOM  (the  synthesis  module)  on  0/P  of  (4),  with  a  unity-gain 
Correction  file.  This  produces  the  synthetic  target  file. 

(6)  Output  this  into  Che  chamber  using  'TESTX4',  while  simultaneously  obtaining 
the  1/3  oct.  data  for  the  actual  chamber  noise  as  in  (3)  above. 

(7)  Obtain  error  (difference)  file  from  (3)  to  (6)  by  DATS/ARITH  and  then  run 
/DIFFER  to  get  a  1st  Correction  file. 

(8)  Use  /XASCOM  om  component-band  files  found  in  (4),  this  time  with  the 
1st  Correction  file.  The  0/P  from  this  operation  is  the  1st  Compensated 
output  file. 

(9)  Repeat  (6)  and  (7)  above  to  examine  the  improvement  in  error  spectrum  and 
to  proceed  to  a  further  iteration  if  necessary  by  repeating  (8),  using  the 
new  (2nd)  Correction  file. 
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Fig  la^b  Third-octave  error  spectra  (jet  aircraft) 
(a)  original  (b)  improved 
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Fig  2  Implementation  of  compensation  orocess 
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Fig  Ai  Periodic  spectra  produced 
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Fig  A4c  Example  of  warped  frequencies  (s-w2  planes) 


Fig  A4d  Warped  frequencies  (s-Wj  planes) 
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Fig  A6  S-plane  iow  pass  to  band-pass  frequency  transformation 


Low  Pass  ({Ia'i/s ec.) 


Fig  AS 


Fig  AO  Final  mapping  of  prototype  low-oass  continuous  filter 
to  digital  bandpass  filter 


